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Summary. — Observing quantum particle creation by black holes (Hawking radi- 
ation) in the astrophysical context is, in ordinary situations, hopeless. Nevertheless 
the Hawking effect, which depends only on kinematical properties of wave prop- 
agation in the presence of horizons, is present also in nongravitational contexts, 
for instance in stationary fluids undergoing supersonic flow. We present results on 
how to observe the analog Hawking radiation in Bose-Einstein condensates by a 
direct measurement of the density correlations due to the phonon pairs (Hawking 
quanta-partner) created by the acoustic horizon. 
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1. Hawking radiation and experimental prospects 

In 1974 Stephen Hawking [1] discovered that due to quantum mechanics black holes 
are not "black" , but emit particles in the form of thermal radiation at a characteristic 
temperature 



where k (= c 3 /4GM for a Schwarzschild black hole) is the horizon's surface gravity. This 
astonishing result was obtained by studying the vacuum of matter fields in the vicinity 
of the event horizon, the black hole's boundary marking the region from where nothing 
(not even light) can escape. 

Quantum vacuum fluctuations are usually pictured in terms of pairs of virtual par- 
ticles being continuously created out of the vacuum, but which almost instantaneously 
(after a time of the order of Planck time tp ~ 1CP 43 s) annihilate and disappear leaving 
no observable effect. The presence of the event horizon can change dramatically this 
picture, since when one member of the pair is created just inside the horizon and the 
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Fig. 1. - Schematic representation of an acoustic black hole in a fluid, vq — Vq(x) and c s = c s (x) 
are, respectively, the fluid and sound velocities. The left region has supersonic flow and the 
right one is subsonic; they are connected through an acoustic horizon. 



other just outside, the one in the interior (the partner) gets trapped and may leave the 
other free to propagate far away where it is detected as a real particle (Hawking quanta). 
This heuristic argument was given by Hawking in [2] to explain the basic mechanism 
responsible for the black hole's emission. 

Because of its physical consequences (namely, the information loss problem [3]) the 
Hawking effect is nowadays considered as a milestone for the construction of a quantum 
theory of gravity. Nevertheless, it is practically impossible to observe it. Since we 
cannot, for obvious reasons, check at the same time the existence of the Hawking quanta 
and of their partners inside the horizon, an experimental verification of the Hawking 
effect requires the detection of the emitted thermal flux at the temperature (1) far from 
the black hole. The bad news is that, unfortunately, in ordinary situations where the 
black holes are created from gravitational collapse of massive stars the Hawking flux at 
Th ~ 10~ 7 K M ^j n < 10 _7 K gets completely overwhelmed by the cosmic microwave 
background radiation at Tqmb ~ 3 K.(}) 

2. — A way out: analog Hawking radiation from acoustic black holes in fluids 

To remedy this disappointing situation, in 1981 Bill Unruh [6] explored an alternative 
route. It was known that the propagation of sound waves in inhomogeneous eulerian fluids 
and that of (massless) scalar fields in curved backgrounds are mathematically equivalent. 
Moreover, sound waves get trapped in regions of supersonic flow in the same way light 
is trapped inside a gravitational black hole. Therefore a stationary fluid configuration 
consisting in a region of subsonic flow and one of supersonic flow connected through an 
acoustic horizon represents an acoustic analog of a gravitational black hole, see Fig. 1. 

Unruh pushed this analogy further at the quantum level. Quantized sound waves 
(phonons) in these acoustic black holes behave like quantized scalar fields in gravitational 



( ) Alternative possibilities to detect it are based on mini black holes (with a much higher Th) 
production, by density fluctuatons in the early Universe [4] and in particle accelerators due to 
the existence of large extra dimensions [5] . 
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black holes. The same analysis Hawking used can then be applied here leading to the 
prediction that acoustic black holes will produce a thermal flux of phonons from their 
acoustic horizon at a temperature given by (for a ID fluid for simplicity) 



(2) 
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where c s is the sound speed, vq the fluid velocity and the derivative is evaluated at the 
acoustic horizon (where c s = \vq\). 

Although we limited here to the case of fluids, these same considerations apply also 
to other condensed matter systems in the long wavelength (hydrodynamic) approxima- 
tion [7]. Besides offering new unexpected possibilities to test the Hawking effect in a 
nongravitational domain, this application also allows to address from first principles a 
delicate point in Hawking's derivation called the transplanckian problem. 

It is well known that in the propagation from the horizon to infinity light signals suffer 
an infinite exponential redshift. Applied to Hawking's analysis, this implies that in order 
for the Hawking quanta to reach the asymptotic region (where they are detected) the 
relevant field modes in the near-horizon region must be characterized by an enormous 
initial frequency, well beyond the Planck scale. This poses a problem because in these 
regimes one should more appropriately use a quantum theory of gravity, the lack of which 
poses serious concerns both on the validity of Hawking's analysis and its final result. 

One can address this problem in the fluid case, since the theory describing scales below 
the interatomic distance (the analogous of Planck scale in gravity) is in many cases well 
known. 

Many directions are being explored, the most recent ones ranging from water tanks 
experiments (stimulated Hawking radiation from white hole flows [8]) to laser pulse 
filaments [9] ( 2 ) and to Bose-Einstein condensates (BECs), the focus of this contribution. 



3. The gravitational analogy in BECs 

Bose-Einstein condensates are ultracold bosonic systems characterized by the fact 
that all constituents share the same quantum state. Perfect condensates exist only at 
absolute zero temperature, while in practice a small non condensed fraction is always 
present. To describe them one writes down the Bose field operator \t = ^>q(1 + 4>), where 
the classical wavefunction v^o describes the macroscopic condensate and the quantum 
operator tp the small noncondensed fraction. In the dilute gas approximation (see for 
instance [11]) satisfies the Gross-Pitaevski equation 

-4^ + Kxt+sl*o| 2 *o , 
2m j 

where m is the constituents mass, V ext the trapping potential and g the nonlinear 
atom-atom interaction constant, while <f> treated as a linear perturbation follows the 



( 2 ) The interpetation of these results is still controversial, see for instance [10]. 
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Bogoliubov-de Gennes equation 

(4) m* =-(^ + * + mc?M + &) 



2m m ^ 



with c s = \p^- the speed of sound. Eq. (3) gives the background while (4) describes the 
quantum fluctuations and is valid at all scales. 

To make contact with the gravitational analogy we need to consider the hydrody- 
namical approximation of the above equations, and this is more easily achieved in the 
density-phase representation in which we write ^ = \fh~e %e . Again we consider the split- 
ting h = no + hi and 9 = 9q + 9\ with no,#o the background and h\,9\ the linear 
fluctuations. Provided the background varies on scales bigger than the healing length 
£ = -5— (the microscopic scale, analogous to the Planck length in gravity), cq. (3) for 
no, 9q reduces to the continuity and Euler equations for irrotational inviscid fluids (the 
analogous of Einstein's field equations) . Similar arguments applied to (4) give two oper- 

atorial equations for the density and phase fluctuations hi = no(0 + (ft), #i = %f in 
the form 



(5) hi = -2 

met 



v d x 9i + d t 9i 



where vq = hd ^° is the condensate velocity (we consider a quasi-fD condensate for 
simplicity), while the equation for 9i decouples and is mathematically equivalent to the 
Klein-Gordon equation for a massless, minimally coupled scalar field 

(6) -^d^^g^dA) = 



in the background metric 

(7) ds 2 = (-(fii -vhdt 2 +2v dtdx + dx 2 + dy 2 + dz 2 ) . 

mc s 

Eq. (6) is the equation Hawking used to quantize a scalar field in the Schwarzschild 
background. The acoustic metric (7) possesses an acoustic horizon located at \vq\ = c s 

with surface gravity n = 2~~ ^"dx"^ \hor, an d unlike its gravitational analogue it does 
not need to have the interior (supersonic) region to terminate into a (nasty) spacclike 
singularity. Since Hawking's analysis only relies on the kinematical properties of the 
horizon, event or acoustic, the prediction that in the formation of an acoustic black hole 
a thermal flux of phonons at the temperature (2) will be produced is then straightforward. 

The experimental realization of an acoustic black hole in a BEC was carried out 
for the first time in in [12]. Because of their extremely low temperatures (of the order 
of Tc ~ lOOnK, much below Tqmb) they offer much more favourable experimental 
conditions for the detection of the Hawking effect with respect to gravity. Indeed, acoustic 
black hole configurations can be constructed for which ~ lOnK, only one order of 
magnitude below Tc- Nevertheless, the signal is tiny and competing effects (the bigger 
background temperature, quantum noise, etc..) still easily overcome the Hawking flux of 
the created phonons making it practically unobservable. 
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4. One step forward: the Hawking effect in density correlations in BECs 

Another experimental advantage of having to do with acoustic black holes instead of 
gravitational ones is that acoustic horizons, unlike event horizons, allow measurements 
to be performed on both sides, the exterior (subsonic) and the interior (supersonic) 
regions. This is very important because we can then try to characterize experimentally 
the Hawking effect through the pair production process responsible for it. This was done 
in [13], to which we refer for more details. 

In [13] the phase fluctuation 6\ was treated as a 2D conformal scalar field propagating 
in the (t, x) section of the acoustic metric (7), for which the two-point function 
is known analytically ( 3 ). This powerful result allowed us to construct in closed form, 
using (5), the equal time normalized density-density correlation function in the Unruh 
state (the state appropriate to describe black hole evaporation). Focussing in particu- 
lar on points situated on either side with respect to the acoustic horizon, we got the 
characteristic stationary behaviour 

(8) O&-4<fc»0~«^-[=<;^ 

Here v < 0, i.e. the flow is from right to left while x = is the acoustic horizon dividing 
the right (x < 0) subsonic and left (x > 0) supersonic regions. The formation of the 
horizon triggers the spontaneous creation of pairs of phonons (of equal and opposite 
frequencies, as required by energy conservation) on both sides of the horizon. In the lab 
they propagate along the lines x = (v + ci)t (< 0) (partner), x' = (v + c r )t (Hawking 
quanta), see Fig. 2. ( 4 ) This happens for all t and leads to a stationary signal correlating 
points x, x' such that 

(9) — ^— = . 

V + Cl V + C r 

The correlator (8) has a peak exactly along this line: this signal is a clear manifestation 
of the pair-production process at the base of the Hawking effect. Quantitative estimates 
based on existing experiments suggest the the magnitude of this peak, of the order of 
10~ 3 , is small but not negligible. 

One can object that this result and the thermal flux prediction, being based on the 
gravitational analogy, are subject to a 'transplanckian' problem analogous to that in 
gravity. Indeed, sound waves propagating from the near horizon region get exponentially 
redshifted, and this implies that the physics of the Hawking effect is in principle sensitive 
to microscopic scales, in particular those smaller than £. In BECs we can study the 
'transplanckian' problem by going back to the basic equations (3) and (4), valid at all 
scales (we do not have the analogous equations in gravity). 

The result (8) served both as a motivation and a guide for a subsequent numerical 
analysis performed within the full theory using Montecarlo techniques [14]. The results, 
shown in Fig. 3, confirmed the existence of the Hawking peak in the microscopic theory. 
A quantitative comparison of the height and width of the peak with (8) showed excellent 



( 3 ) This approximation only slightly overestimates the flux of the emitted phonons. 

( 4 ) Although we are considering its hydrodynamic/relativisistic approximation, the system we 
started with is nonrelavistic and velocities add up following Galileo's law. 
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Fig. 2. - Pair-creation (Hawking quanta - partner) from the black hole horizon. 



agreement with the hydrodynamical prediction in the regime k -C | ■ Outside this regime 
(which in the gravity context would require knowledge of a quantum theory of gravity) 
acoustic black holes still emit at an approximately thermal rate (up to a certain cutoff 
frequency w max ~ j [15]) and the location of the Hawking peak is still given by (9) (for 
analytical results in this regime, see [16] and references therein). 

A nice way to explain why there is no transplanckian problem in BECs is depicted 
in Fig. 4. By tracing backwards in time the trajectories of the Hawking quanta and 
partners, once the near horizon region is reached microscopic physics effects change 
the phonon's dispersion relation which becomes 'superluminal'. This causes the two 
trajectories to stop hovering close to the horizon (where, instead, they would remain 
forever in the hydrodynamical/relativistic approximation) and enter the interior region. 
The resulting blueshift along the trajectories is now finite, of the order of | (we recover 
the infinite blueshift in the hydrodynamcal £ — > limit). 

Another important result of the numerical analysis in [14] is that the Hawking peak of 
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Fig. 3. - Numerical results for the density correlator confirming the existence of the Hawking 
peak (red tongue perpendicular to the main diagonal) in the microscopic BEC theory. 
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Fig. 4. - Because of short-distance effects (scales smaller than f) Hawking quanta and partners 
appear to be produced in the interior of the acoustic black hole, thus resolving the transplanckian 
problem in BECs. 



Fig. 3 is still visible even when a nonzero thermal background (with a temperature bigger 
than Tjj) is present. Thus, unlike the emitted thermal flux (which is easily masked by 
the background ) the correlator signal (8) is the looked for 'smoking gun' of the Hawking 
effect in the experimental search. 

5. Acoustic white holes in BECs: another experimental signature of the 
Hawking effect 

White holes, time reversal of black holes, have not captured much interest in gravity. 
They expel to the exterior everything that lies inside their horizon, exactly the opposite 
behaviour of black holes. They 'emerge' from an initial singularity, a kind of big-bang 
singularity, which prevents physical predictions about them to be made until we'll find 
a way to deal with singularities in the quantum theory. 

This problem is not present in the fluid case, and acoustic black and white holes can 
be studied on the same footing. A schematic representation of an acoustic white hole is 
given in Fig. 5, and is simply obtained by reversing the sign of the fluid velocity Vo in 
Fig. 1 (and also in the acoustic metric (7)). An acoustic white hole is represented by a 
supersonic fluid decelerating until it becomes subsonic, whereas the acoustic black hole is 
the other way around, namely an accelerating fluid turning from subonic to supersonic. 
In both cases, the acoustic horizon lies where the fluid velocity equals the sound speed. 

Spontaneous pair production in stationary backgrounds requires the existence of neg- 
ative energy states. These exist in supersonic regions of both acoustic black holes and 
white holes, indicating that an analog Hawking effect exists in both configurations and 
not just in black holes ( 5 ). Its features in white holes are however totally different. 
Because of the time-reversal property with respect to the black hole case (Fig. 2), the 
trajectories in Fig. 6 show that Hawking quanta and partners now accumulate along the 



( 5 ) Perhaps surprisingly, such an effect exists also for an everywhere supersonic fluid with no 
acoustic horizon, see [16]. What seems to depend on the existence of an acoustic horizon is the 
thermal character of the created flux. 
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Fig. 5. - Schematic representation of an acoustic white hole in a fluid, obtained from Fig. 1 by 
reversing the sign of the fluid velocity vo- 



horizon, and they do this with a diverging blucshiftcd frequency. This implies that hy- 
drodynamical white hole horizons (as well as relativistic ones) are unstable: the Hawking 
effect in white holes crucially depends on the details of the short-distance physics. 

The 'superluminal' modification of phonon's dispersion relation in BECs causes (as 
in the black hole case) departure from the hydrodynamical picture in Fig. 6 when a 
maximum frequency ~ | is reached, allowing both trajectories to enter the horizon (Fig. 

7). 

These arguments indicate that the main Hawking signal in BEC acoustic white holes 
takes place inside the horizon. Analytical results derived within the linear Bogoliubov 
theory show that for x,x' < 0, i.e. inside the horizon, the density-density correlator has 
the typical oscillating behaviour [17] 

(10) G^ H {x,x r ) ~ [Aco$,po{x + x') + Bsinpo(x + x') + Ccospo(x — x')} I e , 

with pa the nontrivial zero mode (~ |) of the 'superluminal' Bogoliubov dispersion 
(present only in the supersonic region), and the overall amplitude I e = f —, formally 
infrared divergent, is regulated by introducing the low- frequency cutoff e = j where 
t measures time elapsed from horizon formation. A closer inspection shows [18] that 
the two-point function in the relevant low-frequency regime factorizes as the product of 
the same function </>o evaluated at the points x and x' ', typical of an emerging classical 
behaviour. It describes the emission of a zero-frequency undulating wave away from the 
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Fig. 6. - Trajectories of Hawking quanta and partners in acoustic white holes. 
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Fig. 7. - Short-distance physics (scales smaller than £) cuase Hawking quanta and partner to 
enter the white hole horizon instead of accumulating along the horizon on both sides (Fig. 6). 



horizon with a macroscopically growing amplitude (i.e. the — term in J e ) fixed by the 
low-frequency spectrum of the spontaneously produced Hawking phonons. It is worth 
to mention that undulations were observed in water tanks experiments in [8], but their 
connection with the Hawking effect was not pointed out. 

6. Conclusions and future prospects 

We have described how the analog Hawking radiation can be measured in Bose- 
Einstein condensates by a direct measurement of the density correlations between the 
spontaneously produced pairs of phonons (Hawking quanta and partners) on both sides 
of the acoustic horizon (such a measurement would be impossible to perform in gravity) . 
This provides the looked for 'smoking gun' of the Hawking effect and makes these type 
of measurements very promising for an experimental verification of the Hawking effect 
in the near future. Methods to amplify this signal have been proposed [19] by letting the 
condensate to expand freely after the acoustic black hole is formed. 

Finally, it is worth to mention that density correlation measurements between atoms 
velocities were performed recently in [20] to observe the creation of correlated excitations 
with equal and opposite momenta in homogeneous condensates undergoing a sudden 
variation of their trapping potential. While this effect can be interpreted as the analog 
of cosmological particle creation, this technique will prove useful to test the Hawking 
effect in stationary supersonic flows. 



This work is part of a research developed within the Project "Acoustic black holcs ; 
of Centro Fermi, in collaboration with Roberto Balbinot. 
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